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Abstract 

A particle rotor model (PRM) with a quasi-proton and a quasi-neutron coupled with a triaxial 
rotor is developed and applied to study chiral doublet bands with configurations of a /in/2 proton 
and a quasi-neutron. With pairing treated by the BCS approximation, the present quasi- 

particle PRM is aimed at simulating one proton and many neutron holes coupled with a triaxial 
rotor. After a detailed analysis of the angular momentum orientations, energy separation between 
the partner bands, and behavior of electromagnetic transitions, for the first time we find aplanar 
rotation or equivalently chiral geometry beyond the usual one proton and one neutron hole coupled 
with a triaxial rotor. 
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I. INTRODUCTION 



Since the pioneering work of Frauendorf and Meng [l| , the phenomenon of chiral rotation 
in atomic nuclei has attracted significant attention. Chirality in nuclei offers direct evidence 
for the existence of stable triaxial nuclear shapes, in which there are a few high-j valence 
particles and a few high-j valence holes. For a triaxially deformed rotational nucleus, the 
collective angular momentum favors alignment along the intermediate axis, which in this 
case has the largest moment of inertia, while the angular momentum vectors of the valence 
particles (holes) favor alignment along the nuclear short (long) axis. The three mutually 
perpendicular angular momenta can be arranged to form two systems with opposite chirality, 
namely left- and right-handedness. They are transformed into each other by the chiral 
operator which combines time reversal and spatial rotation of 180°, \ = T1Z(it). The 
spontaneous breaking of chiral symmetry thus happens in the body-fixed reference frame. 
In the laboratory reference frame, with the restoration of chiral symmetry due to quantum 
tunneling, the so-called chiral doublet bands, a pair of separated AI = 1 bands (normally 
regarded as nearly degenerate) with the same parity, are expected to be observed in triaxial 
nuclei. 



Originally a pair of AI = 1 bands found in Pr with the 7r/in/2(8>z^ii/2 configuration 



a, 



has been reinterpreted in Ref. candidate for chiral doubling. Thereafter, similar 

low-lying doublet bands were reported in 55CS, 57La, and 6iP m N = 75 isotones of 134 Pr, 
and an island of chiral rotation was suggested in the A ~ 130 mass region So far, 
candidate chiral doublet bands have been proposed in a number of odd-odd nuclei in the 



A rsj 130 



a a b, a a b, b , Q, fl Q and A 



A few 



regions k], 14, [lol ]. 

more candidates with more than one valence-particle and hole were also reported in odd 

nnnn n 

A [IS 113, [18|, [19|] and even-even nuclei [20] . 



On the theoretical side, chiral bands were first predicted in the particle-rotor model 
(PRM) and tilted axis cranking (TAC) approach for triaxially deformed nuclei [l|. Numer- 
ous efforts have been devoted to the development of the PRM and TAC approaches. Chiral 
rotation has been studied by the Strutinsky shell correction TAC (SCTAC) method with 
a hybrid potential which combines the spherical Woods-Saxon single-particle energies and 



the deformed part of the Nilsson potential 



21 



221 ] . Recently, chiral TAC solutions have 



also been found in N = 75 isotones within the self-consistent Skyrme Hartree-Fock cranking 
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model 



23 



24( . The cranked relativistic mean field (RMF) theory has been reported only in 



the contexts of principle axis rotation 



25 



261 ] and planar rotation 27J. The generalization 



thereof 
In Ref. 



or searching for chiral solutions, i.e., the aplanar rotation, is still under development. 



2s 



the adiabatic and configuration-fixed constrained triaxial RMF approaches were 
developed to obtain the nuclear potential energy surface with the triaxial degree of freedom 
and the existence of multiple chiral doublets (M%D) was predicted for A ~ 100 mass region 
based on their triaxial deformations and their corresponding proton and neutron configu- 
rations. The advantage of the cranked mean field approach to describe nuclear rotation 
bands is that it can be easily extended to the multi-quasiparticle case. However, the usual 
cranking approach is a semiclassical model, where the total angular momentum is not a good 
quantum number, and the description of quantum tunneling of chiral partners is beyond the 
mean field approximation 28, 3o| . 

In contrast, the PRM is a quantum mechanical model where total angular momentum is 
a good quantum number. The model describes the system in the laboratory reference frame 
and yields directly the energy splitting and tunneling between doublet bands. Chirality for 
nuclei in A ~ 100 and A ~ 130 regions has been studied with the particle-rotor model 



for certain partic 

r 

coupling model [9 



e-hole configurations 
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3jJ , or the core-quasiparticle/core-particle-hole 



321 ] following the Kerman-Klein-Donau-Frauendorf method [33] . Selection 
rules of electromagnetic transitions for chiral doublet bands have been proposed based on a 
simple particle-hole-triaxial rotor model 34]. 

Though various versions of PRM and TAC have been applied to study chiral bands, the 
essential starting point for understanding their properties is based on one particle and one 
hole coupled with a rigid triaxial rotor. Based on this scenario, a set of observable signatures 



35 



36 



371 ] . Critical analyses for 



have been suggested as fingerprints of chiral bands 
the representative cases of candidate chiral bands, 134 Pr in A ~ 130 36(, and 104 ' 106 Rh in 
A ~ 100 37J] have been carried out. It has been found that these candidate chiral bands 
in 134 Pr and 104 Rh do not agree with all of those expected for chiral bands, although these 
candidates have been considered as the best examples of chiral rotation in the A ~ 130 and 
A ~100 mass regions (due to their extremely small level discrepancy between the doublet 
bands). Lifetime measurements are essential for extracting the absolute B(M1) and B(E2) 
transition probabilities, which are critical experimental observables in addition to the level 
energies. Indeed, this has stimulated experimental programs aimed at identifying chiral 
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doublet bands 



38 



39]. 



On the other hand, one should bear in mind that these fingerprints of chiral bands are 
obtained mostly by assuming one proton (neutron) particle and one neutron (proton) hole 
sitting in a high-j shell coupled with a triaxial rotor with 7 = 30°. In a realistic nucleus, it is 
more natural that there will be more than one nucleon in a high j-shell, e.g., the candidate 
chiral doublet bands reported for N = 75 isotones with Z = 55 ( 130 Cs), 57 ( 132 La), 59 
( 134 Pr), 61 ( 136 Pm) and 63 ( 138 Eu), and for Z = 55 (Cs) isotopes with N = 69, 71, 73, 75, 
and 77. The Fermi energy of a proton (neutron) will undoubtedly change with Z (N) in 
these isotones (isotopes). Therefore it is interesting and necessary to investigate the doublet 
bands with valence nucleons sitting in the middle of a high j-shell, or alternatively multi- 
particles sitting in a high j-shell. It is also important to investigate the properties of doublet 
bands as functions of the triaxial deformation degree of freedom. 

To address these issues, in this paper a particle rotor model with a quasi-proton and a 
quasi-neutron coupled with a triaxial rotor is developed and applied to study chiral doublet 
bands with configurations of a hu/2 proton and a hn/2 quasi-neutron. With the pairing cor- 
relations taken into account by the BCS approximation, the configuration of multi-particles 
sitting in a high j-shell can be simulated by adjusting the neutron Fermi energy. Note that 
in a former paper 40j, the present model has been applied to the doublet bands of 126 Cs, 
and good agreement with the data available was obtained, which supports the chiral inter- 
pretation of these doublet bands. Here the formalism is given in detail and the properties 
of doublet bands calculated are presented. The model is introduced in Sec. II. The prop- 
erties of doublet bands thus obtained, such as energy spectra, electromagnetic transitions, 
and the orientation of angular momenta, are discussed in Sec. III. Finally, a summary and 
conclusion is given in Sec. IV. 



II. FORMALISM 



The particle rotor model 4l| for triaxial de 



30 



tion of the odd- A and odd-odd nuclei 
nucleus can be expressed as, 

H = H co n 



'ormed case has been well used for the descrip- 



42 



43|, 



44 



H? 

"intr 



45]. Its Hamiltonian for an odd-odd 

(1) 



TTTl 

12 intr ) 
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where p and n refer to protons and neutrons, respectively. The collective Hamiltonian takes 
the form 



H, 



R 



coll 



2Ji ^ 



U Jpi Jr 

2Ji 



(2) 



i=l i=l 

where 7j, j pi , j ni respectively denote the angular momentum operators for the core, nu- 
cleus, as well as the valence proton and neutron. The moments of inertia for irrotational 

flow are adopted, i.e., Ji = J sin 2 (7 —%) . 

The intrinsic Hamiltonian for valence nucleons is 



Hint? = -^sp + -ffpair = — A)(a+a„ + d^d-p) — — ^(a+a± + a v a v 



A 



(3) 

i/>0 u>0 

where A denotes the Fermi energy, A the pairing gap parameter, and \v) the time-reversal 
state of \u). The single particle energy e v is obtained by the diagonalization of the Hamil- 



tonian H sp . Similar as in Ref. 



30], for a single- j shell, one has 



H sp = +-C { cos j(j d 3 



L sin 7.. 2 , 2 . 



(4) 



where the plus sign refers to a particle, the minus to a hole, and the coefficient C is propor- 
tional to the quadrupole deformation (3 30, |43j. The single particle states are thus written 

as 

«+|0> = £4^, 4|0) = ^(-l^V-n, (5) 

where f2 is the projection of the single-particle angular momentum j along the 3-axis and 
can be restricted to the values • • • , —7/2, —3/2, +1/2, +5/2, • • • due to time-reversal degen- 



eracy. 



44 



45] 



To obtain the PRM solutions, the total Hamiltonian ([I]) must be diagonalized in a com- 
plete basis space, which couples the rotation of the inert core with the intrinsic wave functions 
of valence nucleons. When pairing correlations are neglected, one can construct the so-called 
strong coupling basis as 



\IMKu p u n ) 



1 21+ 1 

2 V 8tt 2 



21 + 1 



zip f^n 



Dm,k«\0) + (-i) I - K D l M,-K<4 p 4jo) 



for K = ±1, ±3, iS- 
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\lMKv p V n ) 



1 121 +1 



21 + 1 



for K 



1571-2 Z^Z^^P 

0,±2,±4< 



(7) 



The restriction on values of K is due to the fact that the basis states are symmetrized 
under the point group D 2 , which leads to K — Q p — Q n in Eq. ([H]) and K — Q p + Q n in 
Eq. (J7|) being an even integer 44j. The matrix elements of Hamiltonian (|2j) and (jlj) can 
be evaluated in the basis ([6]) and (J7|), and then diagonalization gives eigenenergies and 
eigenstates for the PRM Hamiltonian. For a certain spin J, the dimension of the basis space 
will be (1/4) (2/ + l)(2j p + l)(2j n + 1). 

To include pairing effects in the PRM, one should replace the single particle state a+|0) in 
the basis states and ([7]) with the BCS quasiparticle state c*+|0) to obtain a new expansion 
basis, where |0) is the BCS vacuum state. The quasiparticle operators a+ are given by 

,+ , 

(8) 



where u 2 + t> 2 = 1. In this new basis, the wave functions of PRM Hamiltonian are written 

as 

|/M)= Yl (c^\IMKu p u n )+C^ v JIMKu p V n )), (9) 



K,Vp,V n 



in which v p and v n represent the quasiparticle states a^|6) and a^jO) instead. Fur- 
thermore, single-particle energies e v should be replaced by quasiparticle energies e' v = 
\/ (e u — A) 2 + A 2 . The total Hamiltonian then becomes: 



H = H coll + ^2e' Vp (a+ p a Up + a±a Vp ) + 5^e'^(a£ a Vn + a± n a Vn ] 



(10) 



To construct the matrix of the above Hamiltonian, in comparison with the case excluding 
pairing, each single-particle matrix element needs to be multiplied by a pairing factor u^u u + 



43 



45(| . The occupation factor v u of the state v is given by 



e v -\ 



The reduced electromagnetic transition probabilities are defined as 

|2 



4l| 



B(a\,I^I>)=J2\(l'M'\^\IM)f 

HM 1 



(11) 



(12) 
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where a denotes either E or M for electric and magnetic transitions, respectively, A is the 
rank of transition operator, and ls the electromagnetic transition operator. 

For electric quadrupole (E2) processes, the corresponding transition operator is generally 
taken as 

J({E%ii) = j p e (f)r 2 F 2M (^,0)rfr, (13) 

which is proportional to the electric quadrupole tensor operator Q 2/1 with a factor ^/5/167r . 
The quadrupole moments in the laboratory frame and the intrinsic system {Q 2fl ) are 

connected by the relation 

4, = vllQ'20 + + vlU)Q' 22 . (14) 

For stretched E2 transitions, one has 
B(E2,Ia -> I' a) 



1 cQo 
l07T 



KK' 

SHI7 



E ^l^X <™0|I'tf' > cos 7 + ^ «//i22|/'/i'> + (/Ji2 - 2\I'K')) 



+Term2, (15) 

where Qo is the intrinsic charge quadrupole moment and the "Term2" term is the same as 
the first term but with the replacement [y n — * V n ). 

For Ml transitions, the magnetic dipole transition operator can be written as 

JZ{MI,H) = \f^^^.(9p ~ 9R)jpiA + On - 9R)jn^ (16) 

where g p , g n , gn are respectively the effective gyromagnetic ratios for valence proton, valence 
neutron and the collective core, and denotes the spherical tensor in the laboratory frame. 
The Ml reduce transition probability B (Ml) is expressed as 



B(Ml,Ia -> I' a) 



3 

Air 



E r iK r i'K' y P ) n {<) J^LK) 



tiKK' n> p n'„ 



[{IKl^\I'K'){n' p n' n \f^ p n n ) + (-iy-b-i»(I - Klfi\I'K')(Q' p Q' n \f,\ -Q p - fi n )} 
+Term2 + Term3 + Term4, (17) 

where terms "Term2", "Term3", "Term4" are the same as the first term but with the 
replacement (y n — > V n ), [y' n — > V' n ), and [y n —> V n , v' n — > V' n ), respectively. The operator 



in Eq. ljTTI) is given by 

= f(p)(g P - 9r)j p ^ + f(n)(g n - gn)jn^ (18) 

with f(p) and f(n) the pairing factors uu' + vv' for proton and neutron, and j M the rank-1 
spherical tensor in the body-fixed reference frame. 

III. RESULTS AND DISCUSSION 

A. Single particle states in the single- j model 

For the intrinsic Hamiltonian of valence nucleons, we apply the simple single- j model, 



which is a good approximation for high-j intruder orbitals [4l| . The single particle energy e 
corresponding to the Hamiltonian in Eq.(j4]) with l/in/2 j'-shell and C = 0.3 MeV is plotted 
in the upper panel of Fig. [1] as a function of the 7 deformation of the deformed well. This 
C = 0.3 MeV corresponds to a quadrupole deformation of (3 ~ 0.28 for the l/iu/2 subshell 
in the A ~ 130 mass region. When 7 = 0°, i.e., axial symmetrical case, there are six 
discrete states with good quantum number Q (±1/2, ±3/2, • ■ •, ±11/2). These states are 
indexed by v (u=l, 2, 6), and the corresponding energies are denoted by e v . When 
axial symmetry is broken, Q is not a good quantum number, and each single particle state 
v is then a superposition of eigenstates of (j 2 , J3) as in Eq. (131) , and changes smoothly with 
7. It can be clearly seen that for a /in/2 particle, a lower energy is obtained for 7 = 60°, i.e., 
an oblate shape is preferred, while for a hole a prolate shape is preferred. Particularly, for a 
nucleus with a irhn/ 2 8>vh~y 2 configuration, the sum of single particle energies will be fairly 
7 soft with a minimum around 7 = 30°, and the 7 degree of freedom will play an important 
role. Note that the single particle energies for levels 2 and 5 are nearly 7 independent. 

With pairing taken into account by the BCS calculation, the quasiparticle energy e' with 
A = 1.227 MeV and A = 1 MeV is given in the lower panel of Fig. [H The Fermi energy A 
is very close to £5, which is shown by a dashed line in the upper panel. The label for each 
level follows the corresponding one in the upper panel. Since the Fermi energy is A ~ £5, 
the state £5 is now the lowest quasiparticle state located at ~ 1 MeV due to the pairing 
gap A. Another feature is that the quasiparticle energy e' u becomes more 7 soft than the 
corresponding single particle energy e v due to pairing. 
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B. Energy spectra 



In the present PRM, if X n = s§ and A n = for neutron, and X p = Ej and A p = for 



6 



30, 



3j|. In the 



proton, the model discussed here is equivalent to the model in Refs. 
following calculation, X p = E\ and A p = are fixed for the proton, i.e., a pure /in/2 particle 
proton, while X n for the neutron changes from the bottom to the top of the /in/2 shell. The 
coefficient C = 0.3 MeV, which corresponds to a quadrupole deformation of ~ 0.28 for the 
A ~ 130 mass region, and the moment of inertia is J = 30 MeV -1 . For the electromagnetic 
transition probabilities, the intrinsic charge quadrupole momentum Q takes a value of 3.5 
eb, and the ^-factors g p — Qr = 0.7 and g n — 9r = —0.6 are adopted respectively. 

Firstly we investigate the behavior of doublet bands for a nucleus at the deformation 
7 = 30° in which the best chirality of nuclear rotation is expected It should be noted 
that for an asymmetric configuration vr(yf s ^ / 2® l/ ^ii/2 ) the best chirality occurs at a deformation 



30| 



7 = 27° in Ref. 

The calculated rotational spectra for the yrast and yrare bands 1 with the configuration 
7r/i 11 /2®z//i 11 /2 for C = 0.3 MeV and J = 30 MeV -1 , are plotted in Fig.[2J In the calculations, 
the odd proton is fixed to be a pure /in/2 particle, while the odd neutron is treated as a 
BCS quasiparticle with A = 1 MeV and X n = ex, £2, ■ ■ ■ ,£e, respectively. The 1 = 9 state 
energies of the yrast bands are taken as reference points and are separated by 2.0 MeV for 
display. 

From Fig. [2J the energy difference between the yrare and yrast bands increases from 
X n = e e to s 1 . For X n = s 6 , two nearly degenerate bands can be clearly seen, especially for 
the spin interval 13 < / < 17, and the energy difference between the doublet bands is below 
100 keV. This is the classical case where the chiral concept was proposed jl|. When X n = £5, 
the two bands are nearly degenerate with a constant energy separation of ~ 200 keV for the 
spin interval 11 < / < 15 and a gradually increasing energy separation for higher spin. For 
A„ = 64, and £3, the spectra present very similar behavior: (1) at the low spin region / < 14, 
a slight odd-even staggering with opposite phase can be seen for yrast and yrare bands; (2) 
only at low spins (/ = 9, 11, 12), the energy difference of the yrast and yrare states is smaller 



1 In the paper, the yrast band denotes the rotational band which connects the lowest energies with given 
spins / obtained from the present PRM calculations, while the yrare band correspondingly connects the 
second lowest energies. 
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than 250 keV; (3) for spin I > 14, the energy differences between yrast and yrare states 
increase with 7, e.g., ~ 400 keV at 7 = 14 and ~ 700 keV at 7 = 20. For X n = ei, odd-even 
staggering becomes more obvious and the two bands are separated by an average energy of 
~ 800 keV. The case of A n = e 2 is in between that of £3 and that of E\. 

The calculated energy difference E 2 (I) — Ei(I) between yrare and yrast bands at spins 
7 = 12, 13, • • • , 17 as a function of 7 deformation is plotted in Fig. [31 The left panel displays 
the results for a pure hn/ 2 proton particle and a pure /in/2 neutron hole (A n = Eq, A = 0) 
configuration. A symmetric E 2 (I) — Ei(I) curve about 7 = 30° can be seen, which in 
turn is associated with the symmetries of Hamiltonians with respect to 7 = 30°. If we use 
A = 1 MeV instead for neutrons, the symmetry will not strictly hold any more. In detail, 
the smallest energy difference (< 200 keV) takes place at 7 = 30° for all the shown spins, 
and particularly at spins 7 = 15, 17, very good degeneracy is obtained, namely the energy 
differences are 7.2 keV and 4.1 keV, respectively. The energy difference increases when the 
7 degree deviates from 30°, and presents a parabola-like curve. At 7 = 20° and 40°, the 
E2(I) — Ei(I) varies from 100 to 250 keV, while at 7 = 15° and 45°, the difference is around 
450 keV. 

In the right panel of Fig. [3l the results for a pure /in/2 proton particle and a neutron 
quasiparticle with A n = £5 and A = 1 MeV are shown. The E 2 (I) — E\(I) curves are 
still parabola-like, while their minima change with the spin 7. The tendency is that the 7 
deformation with the minimum energy difference decreases with spin. It is noted that for 
7 G (20°, 30°), a near constant energy difference (~ 200-250 keV) is observed. Also, the 
energy difference between the yrast and yrare bands is quite large (exceeds 450 keV) when 
the nuclear triaxiality is not prominent, i.e., 7 < 15° or 7 > 45°. 

Among the candidate chiral doublet bands observed experimentally, there are cases with 



a degeneracy point, e.g., 134 Pr with nhu/ 2 <8> ^11/2 [36[, Rh with t^Qq/? <8> ^11/2 



13 



14j. 



cases with a near constant energy difference, e.g., 126, 128, 130, 132 Cs jol, [l2j and 106 Rh 
Ref. 0] suggests that the near constant energy difference may come from a deviation of 
the core shape from maximum triaxiality and a less favorable treatment for the valence 
proton and neutron as a particle-hole configuration. Here our calculations show that either 
a deviation of the core shape from maximum triaxiality or a deviation of the Fermi energy 
surface from a particle-hole configuration will hinder the level degeneracy and prefer a near 
constant energy difference. 
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It is also demonstrated that the small energy difference between the doublet bands sug- 
gests a triaxiality (20° < 7 < 40°) for the nucleus, in comparison with a difference of more 
than 450 keV for 7 < 15° or 7 > 45°. 



C. Electromagnetic Properties 

Electromagnetic transition probabilities are critical observables which carry important 
information on the nuclear intrinsic structure. Using a simple model for a special configura- 
tion in triaxial odd-odd nuclei, Koike et al. suggested the selection rules for electromagnetic 
transitions in chiral geometry 34|. The selection rules yield staggering of B(M1)/B(E2) 
and B (Ml) i n /B (Ml) ou t values for the partner band as a function of spin J, where B(Ml)i n 
and B (Ml) out refer to reduced electromagnetic probabilities for intraband and interband 
AI = 1 transitions, respectively Such staggering behavior has been regarded as a finger- 
print for chirality in odd-odd triaxial nuclei, and has been extensively used to support the 
declaration of chiral doublet bands [3]. It is also acknowledged that in ideal chiral doublet 
bands the electromagnetic transition probabilities must be identical or, in practice, very sim- 



ilar [36j. In the following, the electromagnetic transition probabilities will be investigated 
with two quasiparticles coupled with the triaxial rotor model in order to study whether such 
behavior of the electromagnetic transition probabilities will be influenced by variations in 
configurations and triaxial deformation. 

Fig.|4]shows the intraband B(E2) and B(M1) values of yrast and yrare bands for different 
neutron Fermi energies with 7 = 30°. In the left panel, when A n = the intraband 
B(E2) values at spins I < 14 are nearly zero. This is because the yrast and yrare bands 
are displaced in energy for the lower spin region due to less defined chiral geometry with 
insufficient collective rotation, and these bands are mainly connected with Ml transitions. 
Note that the interband B(E2) values from yrare band to yrast band are large in this spin 
region. For spin / > 15, the intraband B(E2) values increase gradually. For X n = £5, the 
behavior of intraband B(E2) is similar to the case X n = Eq, which is small at low spins, then 
increases with spin. When A n = £4, or £ 3 , the intraband B(E2) values of the yrare bands 
have large differences in comparison with those of the yrast bands. In general, the B(E2) 
values of the yrast bands are larger than those of the yrare bands, especially for spin I < 17. 
Their values become close to each other with / > 18, where the collective rotation of the 
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deformed core makes an important contribution to the total spin. For X n = £2, or £1, the 
intraband B(E2) values of the yrast band increase with spin regularly, whereas those of the 
yrare band exhibit many irregular oscillations. 

For the intraband B(M1) in the right panel of Fig. HI the values of B(M1) systematically 
reduce as the neutron Fermi energy surface A n , decreases from £ 6 to E\. When A n = £1, the 
Ml transition almost vanishes because both the valence proton and neutron are particle-like 
and their contribution to the magnetic moment is canceled by similar angular momentum 
orientations and different ^-factor signs. Therefore the rotation bands for A n = E\ are mainly 
connected by E2 transitions, and correspond to the so-called doubly decoupled bands. For 
X n = e§ and £5, the intraband B (Ml) values of yrast and yrare bands are similar to each 
other. It can also be seen that the odd-even staggering of I? (Ml) for 7 = 30° is obvious 
when A n = £6, while not so obvious in other cases. 

Fig. [5] shows the B(M1)/B(E2) ratios of yrast and yrare bands for different A n with 
7 = 30°, whereas the ratios of "yrare bands" for X n = £2 and E\ are not presented due 
to their irregular B(E2) values. It is interesting to note that the B(M1)/B(E2) values in 
partner bands are close to each other for A n = £6, £5, and £4, in particular for higher spins, 
although there are noticeable differences respectively in B(E2) and B(M1) values in Fig. HI 
Next we examine the odd-even staggering of B(M1)/B(E2) ratios. For X n = Eq, staggering 
can be found for / > 16 in the partner bands due to the staggering of I? (Ml) values. For 
X n = £5, a delicate staggering for / > 16 can be also seen. Except for / < 18 in the yrare 
band for A n = £3, there is no staggering behavior of the B(M1)/ B(E2) ratios in the other 
yrast and yrare bands. 

Fig. [6] shows that the B(M1)/ B(E2) values for the yrast and the yrare bands at different 
7, i.e., 7 = 15°, 20°, 25°, 35°, 40°, 45°, with neutron Fermi energy X n = e§ (left panel) and 
An = £5 (right panel), respectively. The results with 7 = 30° have been presented in 
Fig. [5j One finds that: (1) for all 7 degrees, the values of B(M1)/B(E2) for the yrast 
bands are close to those in yrare ones not only for X n = Eq, but also for X n = £5 (3 neutron 
holes approximately); (2) the staggering of B(M1)/B(E2) ratios sensitively depends on the 
deformation 7. 
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D. Orientations of Angular Momenta 



The key to the formation of chiral bands in triaxial nuclei is the existence of aplanar total 
angular momentum. Using wave functions obtained from the particle-rotor model, one can 
calculate the expectation values of angular momenta, (JA, (ji), and (Ri). The expectation 
values for the three components of the angular momenta J, R, and j p , j n are given as, 

Ji = yOf), 

Ri = yjiiii-ji) 2 )- (19) 

In Fig. the average contributions of the three components / 1(1 + 1), i — 1, 2, 3 to the 
total angular momentum, are plotted for the yrast band (left panel) and yrare band (right 
panel) with X n changing from e$ to E\. In the calculations, 7 = 30°, 1-axis refers to the 
intermediate axis with the largest moment of inertia, and the 2-, and 3-axis are respectively 
the short and the long axis with 3i = 3?> = 1/^Ji- In all panels, it can be seen that the 
average contributions from I\ increase globally with the total spin, while contributions from 
the other two directions decrease globally. 

For X n = Eq, around I = 13, the contributions from three directions are comparable for 
both yrast and yrare bands. This corresponds to a typical case of aplanar rotation. In fact, 
the contributions to the total angular momentum from all three directions are not negligible 2 
in the spin interval (9 < / < 20). Therefore the aplanar solution is realized for this spin 
interval and chiral doublet bands are expected. The statement is also true for the case 
A n = ^5; with the exception that the contribution from the 3rd component is a little smaller 
compared with the case X n = Eq. As the Fermi energy surface A n decreases, the contribution 
from the 3rd component becomes smaller, the total angular momentum will mainly lie in the 
1-2 plane, and an aplanar rotation of the nucleus becomes a planar one. For both A n , = £4 
and £3, aplanar solutions can only be expected around I ~ 11. For A„ = Ei and E\, there 

— 2 

exist only planar rotations. In Fig. [TJ there are some fluctuations of Ii /I(I + 1) for A n = E\, 
E2, £3, and e^ due to the strong interactions between different bands. 

2 According to quantum physics, the minimum contribution from one direction to an angular momentum / 
is given by the value {\ + 1) — I 2 ] }/I{I + 1), when the angular momentum is perpendicular to this 
direction. 
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The expectation values Ri, and j P i,j n i have been investigated for A n = e§,--- , £1 as 
functions of the spin /. For simplicity, the cases for A n = Eq, e$, and S\ are shown in Figs. [HJ 
[91and[TU| respectively. 

In Fig. [HI for \ n = Eq, similar as in Refs. , the collective angular momentum, and 
valence-proton and -neutron angular momentum align along the intermediate axis (1-), the 
short axis (2-) and the long axis (3-) respectively. Since these three angular momenta are 
mutually perpendicular, a chiral picture results. In Fig. [9], for A n = £5, the configuration 
is similar to one proton plus three neutron holes in a single /in/2 shell. In this case, the 
orientations of R and j p are similar to Fig. [SJ, while the third component of the angular 
momentum j n is reduced. The total angular momentum I is still aplanar, but its inclination 
angle to the 1-2 plane becomes smaller. As the neutron Fermi energy surface decreases, the 
hole-like odd neutron will switch to a particle-like one, and j n will align from the 3-axis to 
the 2-axis. Then the valence proton and neutron both align to the 2-axis, with the collective 
angular momentum along the 1-axis, and together they give the total angular momentum 
in the 1-2 plane. This is a planar solution as shown in Fig. [101 Noted that in all cases 
the expectation values along the 1-axis for j n and j p increase with / due to the rotational 
alignment of odd particles. 

The average core contribution to the total angular momentum can be seen in the upper 
panels of Figs. [HJ [9], [10J In Fig. [HJ we note that the core contribution for the 14 + state 
in both the yrast and the yrare band (R ~ Q.5h) is comparable with the contributions 
from the valence proton and valence neutron. The latter is consistent with the result in 
Ref. 46J. In Figs. [HJand[9l R\ increases by 8h (from ~ 4 to ~ 12h) as the spin / changes 
from 12h to 20h. This demonstrates that the increase of the total angular momentum is 
mainly due to the collective rotation for / > 12. Therefore the transition probabilities 
B(E2, 1 — > I — 2) corresponding to the collective rotation should be large for I > 14. These 
results are consistent with the B(E2) values discussed in Fig. [U For the lower spin region 
near the bandhead (/ < 12), the increase of / comes mainly from the contributions of the 
valence proton and neutron while the contribution from the core stays the same. For A n = E\ 
in Fig. [TOj, the collective angular momentum at low spin range / < 16 exhibits odd-even 
staggering which is consistent with the energy spectrum in Fig. [TJ 
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IV. CONCLUSION 



A particle rotor model with a quasi-proton and a quasi-neutron coupled with a triaxial 
rotor is developed and applied to study chiral doublet bands with configurations of a h\\ii 
proton and a /in/2 quasi-neutron. With pairing correlations taken into account by the BCS 
method, a proton and many neutron holes coupled with a triaxial rotor can be simulated by 
changing the neutron Fermi level from the top hnn orbit Eq to the lowest one E\. 

The energy spectra, electromagnetic properties, as well as the orientations of the angular 
momenta of the doublet bands have been investigated in detail. The results are summarized 
as follows: 

1. Aplanar rotation exists at least for A n , = Eq and X n = £5 in a certain spin interval. The 
contributions from the three axes are comparable to each other for the partner bands. 
This demonstrates that chiral geometry holds even for the valence nucleons deviating 
from a pure particle-hole configuration. 

2. The near constant energy separation ( ~ 200 keV ) between the partner bands, which 
has been observed in many candidate chiral bands experimentally, has been obtained 
for A n = Eq and A„ = £5 for certain spin and deformation 7 intervals. 

3. Either a deviation of the core shape from 7 = 30° or a deviation of the Fermi energy 
surface from a particle-hole configuration will hinder the level degeneracy and prefer 
a near constant energy difference. 

4. For 15° < 7 < 45°, A n , lies between e 6 and e 5 , the B(M1)/B(E2) values together with 
B(E2) , BiM 1) for the yrast bands are close to those in yrare bands, which may hold 
for all chiral bands. 

5. The odd-even staggering of B(M1)/ B(E2) values is strongly influenced by the defor- 
mation 7 as well as the Fermi surface A, which suggest that the odd-even staggering 
of B(M1)/B(E2) values may not be a general feature for the chiral bands. 

With pairing treated by the BCS approximation, the present quasi-particles PRM is 
aimed at simulating one proton and many neutron holes coupled with a triaxial rotor. 
After a detailed analysis of the angular momentum orientations, energy separation between 
the partner bands, and behavior of electromagnetic transitions, it is demonstrated that 
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aplanar rotation or equivalently chiral geometry, does exist beyond the simple one proton 
and one neutron hole coupled with a triaxial rotor. While simulating multiple valence 
particles here by adjusting the Fermi energy, one may argue that the valence particles 
dumped into the BCS vacuum in the present model can not contribute to the moments of 
inertia. However, as the main focus is the nuclei in the A > 100 mass region, the influence 
by such approximation should not result in a serious problem. Of course, a model with 
multi-proton particles (holes) and multi-neutron holes (particles) coupled explicitly with a 
triaxial rotor is necessary. Future work should also be devoted to replace the present single- j 
shell by a more realistic single particle potential, such as, Nilsson potential. 
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FIG. 1: Upper panel: The single particle energy e with single- j Hamiltonian in Eq.(j4|) (j = 
11/2, C = 0.3 MeV) as a function of 7 deformation. The six degenerate levels are respectively 
indicated by 1, 2, • • • , 6, as well as the corresponding third angular momentum components ±1/2, 
±3/2, • • • , ±11/2 at 7 = 0° (which is good quantum number only for 7 = 0°). The dashed line 
indicates the Fermi energy A, which is used to obtain the quasiparticle energy e' in the lower panel. 
Lower panel: Quasiparticle energy e' for the same parameters as a function of 7 deformation. The 
pairing parameters are A = 1.227 MeV, A = 1 MeV. Each level (1, 2, • • • , 6) corresponds to that 
with the same number in the upper panel. 
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Spin I [h] 



FIG. 2: Calculated rotational spectra for the yrast (solid circles) and yrare (open circles) bands 
for the configuration vr/i 11 / 2 <8> ^^11/2 with C = 0.3 MeV, J = 30 MeV -1 , and 7 = 30°. In the 
calculations, the odd proton is fixed to be a pure /in/2 particle, while the odd neutron is treated 
as a BCS quasiparticle with X n = e\,E2,-" >£6j respectively, and A = 1 MeV. The 1 = 9 state 
energies of the yrast bands, assumed to be MeV, are separated by 2.0 MeV for display. 
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FIG. 3: Calculated energy difference i^CO — between yrare and yrast bands at I = 

12, 13, • • • , 17 as a function of 7 deformation. In the calculations, C = 0.3 MeV, J = 30 MeV -1 
and the odd proton is fixed to be a pure h\xl% particle, while the odd neutron is treated as a BCS 
quasiparticle with A = 1 MeV, and A n = eg (Left panel) , £5 (Right panel) . 



21 




FIG. 4: Calculated B(E2) and B{M\) values for the yrast and yrare bands: the same parameters 
as Fig. 2 are used. 
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FIG. 5: Calculated B(M\) / B(E2) values for the yrast and yrare bands: the same parameters as 
Fig. 2 are used. 
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FIG. 6: Calculated B(M\) / B(E2) values for the yrast and yrare bands: the same parameters as 
Fig. 3 are used; X n = e$ (Left panel), £5 (Right panel). 
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FIG. 7: (Color online). For the yrast and yrare bands, the average contribution of three components 
to the total angular momentum (if) / 1(1+1), i = 1, 2, 3 in the intrinsic frame is plotted as a function 
of spin I: the same parameters as Fig. 2 are used. Open squares: 1-axis, open circles: 2-axis, open 
triangles: 3- axis. 
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FIG. 8: (Color online). For the yrast and yrare bands, the expectation values for the three 
components of the collective, odd-neutron, and odd proton angular momenta — defined by Ri = 
yj (Rj), j P i = J {jpi), arid j n { = \J (j^i), (i=l> 2, 3) — are plotted as functions of spin /: the same 
parameters as Fig. 2 are used, except that A n = Eq. Open squares, open circles and open triangles 
correspond to the 1-axis, 2-axis and 3-axis, respectively. 
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FIG. 9: Same as FigjHl except that X n = £5. 
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FIG. 10: Same as FigJHJ except that A n = e\. 
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